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Abstract 

Chen and Nester proposed four boundary expressions for the quasilocal 
quantities using the covariant Hamiltonian formalism. Based on these four 
expressions, there is a simple generalization that one can consider, so that a 
two parameter set of boundary expressions can be constructed. Using these 
modified expressions, a nice result for gravitational energy-momentum can be 
obtained in holonomic frames. 



1 Introduction 

For any field theory, Chen and Nester PJE] proposed certain boundary expressions for 
the quasilocal quantities. They found four special expressions. Two are unique under 
the property of their three requirements: the well defined requirement, the symplectic 
structure requirement and the covariant requirement. The other two correspond to 
boundary conditions imposed on a physically meaningful but non-covariant set of 
variables. Generalizing the latter, in the present paper we make some modification by 
a simple adjustment of the four boundary expressions. Using these four expressions as 
the basis, one can obtain a two parameter set of boundary expressions. The modified 
two parameter set of expressions for general relativity includes one case with a nice 
"positive energy" result in holonomic frames in small regions. 



2 Chen-Nester's 4 boundary expressions 

Following pQ, begin with a first order Lagrangian density (i.e., a 4-form) 

C = dqAp-A(q,p), (1) 

where q and p are the canonical conjugate form fields, A is a potential, we suppose 
that q is a /-form and let e = (—1)^. The corresponding Hamiltonian 3- form is 
defined as 

H(N) := £ N qAp-i N £. (2) 
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Taking the interior product of the Lagrangian density (JTJ 

in£ = iNdq A p — edq A ijyp — iN-h 

= £ N q A p — ei N q A dp — edq A i N p — i N A — d(i N q A p), (3) 

where the Lie derivative or form component is £n '■= i>Nd + di^. We see that the 
Hamiltonian 3-form (density) can be put in the form 

H(N) = N' i H fl + dB(N), (4) 

where 

N^H^ = i N q Adp + edq Ai N p + i N A, (5) 
and the natural boundary term is 

B(N)=i N qAp. (6) 

This is called the boundary expression, because when one integrates the Hamiltonian 
density over a finite region to get the Hamiltonian, the boundary term leads to an 
integral over the boundary of the region. However this boundary term is not unique, 
one can add to it whatever he likes without changing the Hamilton equations. Indeed 
this boundary term can be removed by introducing the new Hamiltonian 3-form as 
follows 

H'{N) =H{N) + d{-i N qAp) =N»H IX . (7) 
The variation of the Hamiltonian (jSJ) has the form (for details see [3]) 

8H\N) = -i N (F.E.) — 5q A £ N p + £ N q A dp + dC(N), (8) 

where the field equation term 

F.E. :=5qA^f + ^fA5p, (9) 
oq dp 

is proportional to the first order Lagrangian field equations and can be assumed to 
vanish. The boundary variation term is 

C(N) = -i N q A5p + e5q Ai N p. (10) 

This boundary variation term cannot be removed because it comes from 8H directly. 
Boundary conditions are obtained by requiring that the boundary term in the vari- 
ation of the Hamiltonian vanish. We have to add a suitable boundary term to the 
Hamiltonian 3-form (jSJ) 

H'(N) -> H k (N) = N»H„ + dB k (N), (11) 
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to modify the variational boundary term, in order to get nice components like i]y(8qA 
Ap) or i N (AqASp). In order to achieve such forms, there are just four simple boundary 
expressions that can be added, as can be seen by referring to (|10p. The variation of 



the four Hamiltonians including these different boundary expressions are 

5H q (N) = K + di N (5qA Ap) = K + d(i N 5q A Ap + e5q A i N Ap), (12) 

5H P {N) = K-di N (Aq A5p) = K - d(i N Aq A 5p + tAq Ai N 5p), (13) 

5H d (N) = K + d(-i N Aq A 5p + eSq A i N Ap), (14) 

5H C (N) = K + d(i N 5qAAp-tAqAi N 5p), (15) 

where K is defined as 

K := -i N (F.E.) -5qA £ N p+ £ N qA5p. (16) 

Thus we recover the Chen-Nester 4 boundary expressions 

B q (N) = i N q A Ap-tAq Ai N p, (17) 

B P (N) = i N qAAp-eAqAi N p, (18) 

B d {N) = i N q A Ap-tAq Ai N p, (19) 

B C (N) = i N q A Ap - tAq Ai N p, (20) 



(refer to [3 E] for an explanation of the terminology) where B d stands for ^dynamics, 
B c for £> cons t r aint5 here Aq = q — q, Ap = p — p, q and p are the background reference 
values. Note that the above four boundary expressions can be written in a compact 
form 

B kl ,k 2 (N) = B P (N) + hi N Aq A Ap + tk 2 Aq A i N Ap, (21) 
where (hi, k 2 ) = (0, 0), (0, 1), (1, 0) and (1, 1). In detail 

B , = B P , B 0>1 = B d , B 1>0 = B C , B 1A = B q . (22) 

The boundary conditions associated with B q and B p are the simplest. For the case 
of B q , we want the variation boundary term to satisfy iN^q A Ap) = 0. There are 
two ways to achieve this requirement. First control q, then 5q = 0. The second is 
to freely vary q and then Sq can be anything which implies that we want Ap = 0. 
This is the "natural boundary condition" , it forces p = p, where p is the background 
reference. Similarly for the variation boundary term i^(Aq A 5p) = 0. 

A simple example of the idea is the one dimensional spring with a mass block, 
control the external force 5F and the length of the spring Ax changes or response. 
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Likewise, control the position of the mass block by 8x, there must be a force AF from 
outside the system as response. However the variation boundary expressions from B& 
and B c are more complicated, because they both have mixed control-g, p-response 
and control-p, g-response type terms at the same time. 

3 The modification of the Chen-Nester 4 bound- 
ary expressions 

Consider equation (JT2J), as long as the variational boundary terms pieces i^q A Ap 
and 8q A InAp vanish. It seems that there is no restriction on the magnitude. In 
order words, it can be any real number that one can impose. Consider the variation 
of the boundary expression in a simple way 

SH q = K + d(k 1 i N 5qAAp + ek 2 5q A i N Ap) . (23) 

However, one could construct other similar forms proceeding from 57i p , 57id or 57i c . 
The general form of the variational Hamiltonian of the four different expressions, the 
modified equations (fT2"|) to (fTB^) in one general form all together is 

5H = K + hi N 5q A Ap + b 2 5q A i N Ap + b 3 i N Ap A 5q + b A Ap A i N 5q, (24) 

where hi to 64 are real numbers (it will turn out that they cannot be all independent). 
All of these extra terms can be obtained by adding a suitable multiple of i^Ag A 
Ap or Ag A i^Ap. Actually one can rewrite the modified Chen-Nester 4 boundary 
expressions in a simple way by adding ciixAq AAp + ec 2 Aq AirjAp for the individual 
boundary expressions, where Ci,c 2 G 3?. In detail, the 4 equivalent forms are 

B q = ^Nq A Ap — eAq AtNP + Cii^Aq A Ap + ec 2 Ag A iNAp, (25) 

Bp = ixq A Ap — eAq AinP + C\iNAq A Ap + ec 2 Aq AinAp, (26) 

£>d = iNq A Ap — eAq AinP + CiixAq A Ap + ec 2 Aq AixAp, (27) 

&c = ^nQ A Ap — eAq A i^p + Cii^Aq A Ap + ec 2 Ag A i^Ap. (28) 

They are equivalent because a suitable choice for c\, c 2 in (|23j) can reproduce (|2l)j) to 
(pZSj) with new values for c\, c 2 . For example in (|23j). when C\ — > (ci — 1) gives (|2*7j) . 
In short, one can rewrite the above four expressions in a compact form 

B CliC2 (N) = B P (N) + c x i N Aq A Ap + ec 2 Ag A i N Ap. (29) 
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If the modification can be adjusted by the above 2 parameters of c% and c 2 , then what 
is the role they are playing? Let's consider the asymptotic order magnitude of the 
variables q, q, p and p as follows 

9«0(1)|0Q, 3«0(1), P~°(^)i (3°) 

Since 

AgAAp^o(^), (31) 

then the terms with the coefficients of C\ and c 2 asymptotically have 0(l/r 3 ) fall 
off. Consequently they do not change the dominate main asymptotic value. This 
shows that the modified expression of Cii^Aq A Ap + ec 2 Aq Ai^Ap adjusts the higher 
order terms; but not the lowest order terms, the dominate contribution. Indeed the 
modified expressions are really doing a modified job. Consider the variation of the 
Hamiltonians 

5H q = K + d[(ci + l)i N 5q A Ap + e(c 2 + l)Sq A i N Ap + c x i N Aq A 8p + ec 2 Aq A i N Sp], (32) 
S7i p = K + d[ciiN$q A Ap, +ec 2 5q A iNAp + (ci — l^Ag A 5p + e(c 2 — l)Aq A i^Sp], (33) 
5H d = K + d[cii N 5qAAp + t(c 2 + l)5q Ai N Ap+ (c x - l)i N Aq A 5p + ec 2 Aq A i N 5p],(3A) 
5H C — K + d[(ci + l)iAr5g A Ap + ec 2 5q A i N Ap + C\i N Aq A 5p + e(c 2 - l)Aq A i N 8p].(35) 

Basically all of these four boundary expressions both have terms of the form control-g, 
p-response and control-p, g-response. The boundary conditions associated with these 
four expressions are somewhat like the expressions of £>d and B c . 



4 On the physical interpretation of the modified 
Chen-Nester's 4 boundary expressions 

The modified expressions look a little strange, is this just a mathematical complete- 
ness game or they do have some real physical contribution? What is the idea here? 
Does it have any application? In the Sturm-Liouville theory of 2nd order linear 
differential equation. The necessary condition for self-adjoint is (see e.g. j3j) 

p(yiy' 2 -y2y[)\ b a = o, (36) 

if p(a) ^ ^ p(b) one needs to impose boundary conditions on the solutions. The 
simplest is Dirichlet y(a) = = y(b) or Neumann y'(a) = = y'(b). But a more 
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general choice is to take some linear combination ay + f3y' to vanish on the boundary. 
This is an example of the sort of thing we are considering. Such a boundary condition 
is appropriate for some practical applications such as an elastic cord which is held 
and shook. Then the slope at the end point will be proportional to the end point 
amplitude. 

But when we consider field theory, the idea seems not very suitable. For exam- 
ple for an electrostatic system, such as a parallel plate capacitor, we know how to 
physically impose Dirichlet or Neumann boundary conditions, by fixing the voltage 
or the charge density, but it seems not at all simple to construct a device to fix some 
linear combination of the voltage and charge density, and it is hard to imagine why 
one would want to so such a thing. 

Turning from the well understood electrodynamics boundary problem to the not 
yet so well understood gravitational field boundary value problem, we might at first 
think that the new combined boundary condition expressions would have little use. 
However a little more thought suggests a different view. Unlike the electrodynamic 
case we do not know physically how to arrange for Dirichlet or Neumann boundary 
conditions. When I imagine trying to change the metric on the boundary by moving 
a mass outside, it seems that there will be associated changes in the normal derivative 
of the metric. Maybe it is physically easier to fix a combination of the metric and its 
derivatives. In any case, we found an interesting analytic application for our idea to 
gravitational energy. 



5 Gravitation application of the new quasilocal ex- 
pressions 

For an application of the new modified Chen-Nester expressions consider gravity 
theory. With k = 8nG, let 

g^TV P^^Va P - (37) 

Here T is the connection one form and i] al3 = *{6 a A 9@) where 6 a is the coframe. For 
gravity following [2J E] , rewrite (J29J) as 

2kB CuC2 (N) = 2kB p (N) + c^VAr^ A ArjJ - c 2 AY a p A i N Ar]J, (38) 
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where 

2kB p {N) = AT a p A i N riJ + DpJFArjJ. (39) 

Since we are interested in the energy-momentum components, the DN terms in (|39p 
can be ignored (i.e. we can presume DN = 0). Rewriting (|39)l 

2kB p {N) = AY a p A i NV J, (40) 

Using (gOJ), rewrite (JSE) 

2kB C1iC3 (JV) = Ar^ A ^/ + aiNAVp A At// - c 2 Ar a ^ A i N AnJ. (41) 

For all ci and C2, this expression will give good values at the spatial infinity limit. 
We can find preferred values for c\ and c 2 by considering the small region limit. 

Consider first using Riemann normal coordinates and the adapted orthonormal 
frames. As mentioned in [7j, if we use orthonormal frames and the appropriate ref- 
erence (r p = 0) with N a = constant within matter we get the expected material 
energy-momentum tensor. In vacuum to lowest non-vanishing order we get the Bel- 
Robinson tensor with a positive coefficient only for the case c\ = c 2 = 0. As explained 
in more detail in |6 , , we should get Bel- Robinson tensor because it has positive energy. 

On the other hand we can consider the same formal expression, but using holo- 
nomic variables and reference so that r p = 0, N a = constant in the holonomic frame. 
Again to lowest order inside matter we get the desired material limit. In vacuum, 
however, the basic term B P (N) (it is just the quasilocal expression of KBLB [S] in the 
limit it reduces to the Freud superpotential which gives the Einstein pseudotensor) 
gives the value jg(AB al 3\ a — S a {3\ (T )x x x cr as has long been known This is not so 
good because the vacuum energy is not positive. Can some choices of ci, c 2 save the 
day? The answer is yes. Rewrite (J4*U|) 

2kB p (N) = r a p A i N r)J = - l -V^~gN a U a ^e, u , (42) 
where the Freud superpotential is 

^a [H = -V^g p<T V\p8^ a - (43) 

Consider the C\ term in (J41|) 

i N T% A ArjJ = - l -^N a {h^r aiT - hTT^e^ (44) 
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and the c 2 term 

1 



an } 



r a p Ai N Ar } J = ~J=jjN a { -{8^V\ p -5 v Ji^V^ p ) ) e^. (45) 

+(<^rV - <^rV) 



2 

Therefore (|41|) can be rewritten as 



2kB Ci>C2 (N) = -lN a U^% u + l Cl ^N a (h^T\ n -h^ an )e, 



1 



\'<> L cm - 1 - c 



-c 2 V^iV^ -(^r%-^OT%) (46) 
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+(^rV-^rV) 



From now on, the weighting factor sj— g will be dropped for convenience. Using 
the pseudotensor can be obtained as 

( —Tj H + ^f/iW" _/,"r^ "1 1 

~ "\ -c 2 (rp m - /i^r^ - w% + + ^r\" - V) j 

Inside matter at the origin a short calculation gives 

2^/(0) = 2(7/(0) = 2^/(0). (48) 

Just what we expect from the equivalence principle. In detail the energy density 
inside matter at the origin, the zeroth order term is 

e = -f„°(0) = = -To°(0) = P , (49) 

K 

where k = 8irG and p is the mass-energy density. The momentum density is 

V k = -t\ = -9jl = _ T o fc . (50) 

K 

At the origin in vacuum, the zeroth and the first derivative are 

tj(0) = = <V/(0). (51) 

The first non- vanishing contribution appears at 2nd order. The non-vanishing second 
derivatives in vacuum at the origin, after a little lengthy computation, are 

d 2 ^t af3 (ci, c 2 ) = \ {(4 + ci - 5c 2 )B a fr lu - (1 - 2ci + c 2 )S a/3fU/ + (a - 3c 2 )K af3flu } , 

(52) 
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where the Bel- Robinson tensor B a/3fiu , tensors S a p^ v and K a/3fll/ are defined as follows 

B a P[iv '■— RaX^Rf^ V U + RaXuaR^ff ~ T.9af3Q ' iivR\<r prR^ (53) 

O 

S a (3^,u '■= Ra[i\aR/3v Xa + RavXaR/l^ + ^5 'a/35 'fiuRxa prR X ' Tf " ~> (54) 

3 

Ru /i -gapg^RxaprR^ ■ (55) 

o 

The Bel- Robinson tensor has many nice properties [5] including energy positivity [Hj. 
Consider (jH2I), when (ci, c 2 ) = (0, 0), (0, 1), (1, 0) and (1, 1), they are classified as the 
original Chen-Nester four boundary expressions. The results from [7] are 

0VMO,O) = -(IB^-S^), (56) 

9 2 ^tap^O, 1) = —-(Bapnv + ISaPnv + 3K a p^), (57) 

y 

d 2 p U t a p(l, 0) = -{SBap^ + Saf3[j,v + Kap^), (58) 

y 

2 

#V*a,a(l, 1) = —z K °$v»- (59) 
y 

None are of the desired pure Bel-Robinson form. The general form of the Taylor 
expansion for the Chen-Nester four expressions in compact form is 

f (4 + A* -5*2)5^ ] 
fc 2 ) = 2GJ + — { -(1 - 2/d + A; 2 )5/^ I + 0(Ricci, x) + O(x 3 ),(60) 
18 [ +(h - 3k 2 )K a ^ K J 

Consider (J52)l again, we want the coefficients of and K a p^ u to vanish. Taking 

(ci,c 2 ) = (§,|) gives 

,2 . /3 1\ 2 



^ ( g, g J — -B a pnv (61) 

This result is good, because it only contains the Bel-Robinson tensor. Consequently 
the small region energy will be positive. The general form of the Taylor expansion of 
the expression is 

tj (~, I) = 2GJ + -S/^xV + 0(Ricci, x) + 0(x 3 ). (62) 
V5 5/ 5 

The corresponding four momentum within a small coordinate sphere is 

F * (I = l**^* 1 ** = -^g b ^' = -m, B ~» < 63 > 



where B^ 1 = B^ Q0Q , P M = (-E, Pi) and energy E > 0. Alternatively the four 
momentum can be written in terms of the electric and magnetic parts as follows 

P, (| I) = ~^(E ab E ab + H ab H ab , 2e c ab E ad H b d ), (64) 
where -B^ooo = {E a bE ab + H ab H ab , 2e c ab E ad H b d ) , here E ab = R 0a ob, H ab = ^C 0amn e b mn . 

In an earlier work we have constructed a 10 parameter class of new superpoten- 
tials that give rise to pseudotensors which have positive Bel-Robinson small vacuum 
limit [TQ1E]- But they all seemed very artificial. Our c\, C2 expressions are special 
cases. However they, in contrast, have clear meanings in terms of the boundary con- 
ditions in the Hamiltonian formalism. Here we found a simple specific value for c±, 
C2 which produces the desired positive vacuum result. There is one unique holonomic 
Hamiltonian boundary expression with this property. It is not one of the 4 previously 
considered quasilocal expression but rather a certain specific combination which cor- 
responds to fixing a special combination of Dirichlet and Neumann type boundary 
conditions. Obtained from the boundary term in the variation of the Hamiltonian, 
namely 

j (J^O A At?/ - -51^ A i N A V J - ^i N AT a p A 8 V J + t AV a p A i N 5 Va ^ . 

(65) 

6 Conclusion 

Chen and Nester proposed the four boundary expressions for the quasilocal quantities 
by using the covariant Hamiltonian formalism. Two of them are unique and corre- 
spond to imposing boundary conditions on a covariant combination, the other two are 
non- covariant. Generalizing the latter two cases, one can make some modification by 
a simple adjustment of the four boundary expressions. The basis is using their four 
expressions. Applying the modified expressions to gravity, we found using holonomic 
frames a nice Bel- Robinson result in vacuum for certain specific c±, 0%. 
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